Cavity QED characterization of many-body atomic states in double-well potentials — 

The role of dissipation 
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When an incident light beam is scattered off a sample of ultracold atoms trapped in an optical 
lattice, the statistical properties of the retro-reflected field contain information about the quantum 
state of the atoms, and permit for example to distinguish between atoms in a superfluid and a Mott 
insulator state. This paper extends our previous analysis of this problem to include the effects of 
cavity damping. We use a Monte Carlo wave function method to determine the two-time correlation 
function and time-dependent physical spectrum of the retro-reflected field in the case of a simple 
two-well lattice. We also analyze quantitatively the entanglement between the atoms and the light 
field for atoms in a Mott insulator and a superfluid state. 
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I. INTRODUCTION 

Ultracold atoms in optical lattices provide a remark- 
able test system to simulate a number of situations 
in condensed matter physics [l| under exquisitely con- 
trolled conditions. The first example along these lines 
was the superfluid to Mott insulator transition 0, Q in 
bosonic systems, but the list of strongly correlated con- 
densed matter systems that can be simulated by atoms 
or molecules [1, B 13 in optical lattices has continued 
to grow since these pioneering experiments. Examples 
include, but are not limited to the Bose-Hubbard and 
Fermi-Hubbard models [p, spin systems and the 
Anderson lattice model In other examples, rotat- 

ing lattices are expected to lead to the realization of 
analogs of the quantum Hall effect, and random lattices 
have recently been used to study Anderson localization 
in atomic systems jT3| . 

The ability to characterize the many-particle state of 
the atomic fields in optical lattices is of course central 
to the realization of these experiments. Some schemes 
are designed for probing those states nondemolitionly 
[13, [l^. In a recent paper [3, we proposed an opti- 
cal scheme based on the diffraction of a quantized light 
field off the atomic sample to probe the number statis- 
tics of the matter-wave field. The basic idea is to use 
two light fields counterpropagating in a high-Q ring cav- 
ity and coupled via Bragg scattering off the atoms, a 
technique that is analogous to the Bragg reflection of X- 
rays off a crystal, but operating in the quantum regime. 
We found that the dynamics of the light field strongly 
depends on the manybody state of the atomic field as 
well as on the lattice spacing. Specifically, the statistical 
properties of the Bragg-reflected light field for atoms in a 
Mott insulator state and for a superfluid described both 
in terms of a number-conserving state and a mean-field 
coherent state were found to provide a clear signature of 
the state of the atomic field. 

The present paper extends these results to include the 
effects of cavity damping. We use a Monte Carlo wave 



function method [17| to determine the two-time corre- 
lation functions and (time-dependent) physical spectrum 
p^] of the retro-reflected field in the case of a simple two- 
well lattice. We also analyze quantitatively the entangle- 
ment between the atoms and the light field for atoms 
in a Mott insulator and a superfluid state. Our main re- 
sult is that even in the presence of dissipation these quan- 
tities allow to easily distinguish between a Mott insulator 
and a superfluid atomic state, and in addition permit to 
decide between two familiar descriptions of that state. 

Section II describes the main elements of our model, 
introducing in particular the effective non-hermitian 
Hamiltonian required in the Monte Carlo wave function 
simulations of the problem. The intensity reflected by 
the atoms via Bragg scattering is discussed in section 
III, and section IV presents the time-dependent physical 
spectrum of the reflected light. The quantum entangle- 
ment that may develop as a result of Bragg diffraction is 
quantified in section V in terms of the logarithmic nega- 
tivity. Finally, section VI is a conclusion and outlook. 



II. 



MODEL 



We consider a sample of ultracold bosonic two-level 
atoms @ with transition frequency uja trapped in the 
lowest Bloch band of a one-dimensional double-well po- 
tential [13, [m, [l^l placed inside an optical ring resonator. 
The atoms are driven by two counterpropagating cavity 
modes of wave vectors ±fc and frequency ojk = kc, see 
Fig. [TJ This interaction is described by the dipole inter- 
action Hamiltonian 



Vd = y" dxi)l{x){e\d ■ E{x)\g)i)g{x) + h.. 



(1) 



where d is the dipole moment of the transition, 'Ei(x) 
is the electric field operator, and tpe and ipg are field 
operators describing atoms in their excited and ground 
electronic states |e) and \g), respectively. We assume that 
the optical fields are sufficiently detuned from the atomic 
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transition frequency Ua that the excited electronic state 
can be adiabatically ehminated. 

We proceed by expanding ipe and ipg on the Wannier 
basis of the lowest Bloch band as 



(2) 



where Cm^ and Cm' are the annihilation operators for ex- 
cited and ground state atoms in the m*^ well, and ipm 

/,(9) 



and tpm are the corresponding wave functions. Introduc- 

(3) 



ing also the operator 

N{d) = no + fiie 

where 



2ikd 



(4) 



and d is the well separation, the atom-field system is 
easily seen to be described by the Hamiltonian [16| 



H = hujkd\ak + f>'9[N{Q){a\dk + aLfcQ-fc) 



-I- N{d)a^_^ak+ N{-d)ala^k 



(5) 



where 
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2p^ 



dx8±k{x)i^i'^*{x)4'\x) 



and p is the dipole matrix element of the atomic transi- 
tion. We drop the label [g) in the following for notational 
clarity since no confusion is possible once the excited elec- 
tronic state has been eliminated. 

Cavity damping is treated in the usual way by cou- 
pling the cavity modes to a Markovian reservoir of har- 
monic oscillators of frequencies {i^q\, with the interaction 
Hamiltonian 



h{ak + a-fc) ^ gqb\ + h.c. 



(6) 



Finally, the -\-k cavity mode is driven by an oscillating 
classical current of amplitude rj and frequency uj = oJk, 



Vp = h-qe 



(7) 



The numerical studies that are the subject of this paper 
are conveniently carried out using a Monte Carlo wave 
function approach [l^. In this case, the coupling of the 
system to the reservoir is described in terms of the effec- 
tive non-hermitian Hamiltonian 



ih 



(8) 



where 




Figure 1: Atoms trapped in an optical lattice with lattice con- 
stant d and interacting with two counter-propagating modes 
with wave vectors ±fc in a ring resonator. F: filter. D: Detec- 
tor. 



and 7 is the familiar Wigner-Weisskopf decay rate result- 
ing from the coupling of the cavity modes to the Marko- 
vian reservoir. For the specific situation at hand this 
gives, in an interaction picture with respect to the Hamil- 
tonian Hq = ^[^fe + 5^(0)]afeafe; 



H, 



aft 



III. 



ihj a\ak + a\a~k 



REFLECTED INTENSITY 



(10) 



The light intensity transmitted and reflected by the 
trapped atoms depends strongly on the state of the atoms 
as well as on the well separation d. Throughout this pa- 
per we assume that the incident +k mode is initially in 
a coherent state (with amplitude a — \f2 in our simula- 
tions) and the reflected — fc mode is in the vacuum state, 
and we concentrate on the special cases of well separa- 
tions d = A/4 and A/2, for which the Hamiltonian TJeff 
becomes 



+ h^ (^ate*f(""+"i)* + /i.e.) 



(11) 



Here the 'plus' sign corresponds to the d ~ \/2 case and 
the 'minus' sign to the d = A/4 case. 



A. Mott insulator 

When in a Mott insulator state the atomic field has a 
well-defined number of atoms in each well, and thus can 
be described by a product of Fock states. Atomic Fock 
states are eigenstates of the effective Hamiltonian (|lip . 
that is, the operators no and hi are constants of motion 
with eigenvalues no and ni, and it is possible to replace 
the operators by these eigenvalues in Eq. pT|) . 

It is known [23| that if a system of coupled harmonic 
oscillators satisfies Heisenberg equations of motion that 
can be expressed as 



27 a.k, 



C2 



27 a_fc 



(9) 



'a,=F,{{ak{i)},t), j = l 



(12) 
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Figure 2: Single trajectory (or average reflected intensities 
{n-k{t))) for atoms in an Mott insulator state with 3 atoms 
per well for (a) d = A/4 and (b) d = A/2, rj = 1.5g and 
7 = 0.95. 



Figure 3: Reflected intensity {n_fc(t)) for atoms initially in 
the superfluid state Ii/jspi) with 4 atoms on average per well 
and d — A/4, (a) typical single trajectory, (b) average over 
100 trajectories, for 77 = 1.5g and 7 = 0.9g. 



where the functions Fj may depend explicitly on time, 
then if the oscillators are initially in a coherent state 
they will remain in a coherent state for all times. This is 
the case for the situation at hand. Hence the two modes 
of the light field are in coherent states whose amplitudes 
exhibit damped oscillations due to the combined effects 
of photon exchange between the incident and reflecting 
modes and of cavity decay. As is to be expected, the 
oscillation frequency is proportional to the total number 
of atoms for the d — X/2 case and to the difference in 
the populations of the two wells for d = A/4. The dif- 
ference between the oscillation frequencies in these two 
cases permits therefore a full determination of the well 
populations. 

Figure [5] illustrates the time dependence of the re- 
flected intensity for equal well populations and for the 
two special well separations. Different atomic popula- 
tions in the two wells also result in an oscillating signal 
for the "destructive well separation" d = A/4, but they 
are is much weaker and slower than that in the "con- 
structive case" d = A/2. 

We mentioned that the light field remains in a coherent 
state at all times. These are eigenstates of the annihila- 
tion operator, hence quantum jumps do not affect them, 
and all Monte Carlo Wave function trajectories are iden- 
tical in this case. 



B. Superfluid 

There are two particularly simple ways to approximate 
the state of the atoms in a superfluid state. The first one 
is essentially a mean-field approach that assumes that the 
atoms in the two wells are in coherent states, 

IV'SFi) = |ao,ai) = X! ^no,ni\nQ,ni), (13) 



whereas the second description accounts for the fixed to- 
tal number TV of atoms and describes their state as 

\^SF2) = AA-i( 4+ct)^|0,0) 

N 

= Y,bno\no,N-no). (14) 

no 

This section compares the reflected optical field corre- 
sponding to these two descriptions, assuming as before 
that the incident light field is initially in a coherent state 
and the reflected field is in a vacuum state. 

Consider first the coherent state description of 
Eq. (flS]) . From the discussion of the Mott insulator case, 
we know that for fixed atom numbers in the two wells 
the state of the field remains a coherent state, 

la, 0; no, Til) \ak{no,ni,d,t),a-kina,ni,d,t);nQ,ni), 

so that without quantum jumps between times and t 

\m) = 

Cno.ni \ak{no, ni,d, t),a^k{no, ni,d, t); no, ni). 

(15) 

It is apparent from this expression that the atoms and the 
light field become entangled, a point to which we return 
in section V. Also, since the optical field does not gener- 
ally remain in a coherent state in this case, the quantum 
jumps resulting from the loss of photons due to dissipa- 
tion become apparent in the single trajectories, see Fig.[3l 
Comparing the present case to the situation for a Mott 
insulator state, we observe that there are no contribu- 
tions to the backscattered intensity from situations with 
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Figure 4: Reflected intensity (n^kif)) for atoms initiaUy in the 
superfluid state |V'sf2) with a total of 8 atoms and d — A/4, 
(a) typical single trajectory; (b) average over 100 trajectories, 
for 77 = 1.5(7 a-nd 7 — 0.9g. 

equal populations in the two wells, since the coupling 
between the incident and reflected fields is zero in these 
cases. Also, the oscillatory character of the back-reflected 
light is now largely washed out by the uncertainty of the 
atom numbers in the two wells. 

The number-conserving description p4p of the super- 
fluid state leads to similar results, although the details 
of the dynamics are slightly different because of the fixed 
total number of atoms which results in less atom num- 
ber uncertainty, see Fig. [H In practice, it is therefore 
not expected that the backscattered intensity will allow 
to unambiguously distinguish between the two descrip- 
tions. 



Figure 5: Two-time correlation function G{t+T, t) in the Mott 
insulator regime, with 6 atoms in well and 2 in well 1. Well 
separation is d = A/4, the pump constant is ?7 = O.lg and the 
decay constant is 7 = 0.5(?. 




Figure 6: Two-time correlation function G{t -\- r, t) for atoms 
initially in the superfluid state IV'sfi) with a mean atom num- 
ber of 4 per well and d = A/4, 77 = O.lp and 7 = O.Sg. The 
curve is the average over 50 trajectories at time t and 3 tra- 
jectories at time r. 



IV. PHYSICAL SPECTRUM 

We now turn to the analysis of the physical spectrum of 
the reflected optical field, which also provides a signature 
of the state of the atomic field. 

Because the process under study is not stationary, 
the usual spectrum obtained from the Wiener-Khintchine 
theorem is not appropriate, and we use instead the time- 
dependent physical spectrum of Ref . |18i] . We recall that 
this spectrum is defined as 

Jo Jo 

X (al^ih)d-kit2))dhdt2, (16) 

where H{t,uj,T) is the response function of the filter, uj 
its setting frequency, and T its bandwidth, see Fig. [1] 
Following Ref. [l^ i we assume that it is constant for the 
frequency range of interest, and choose the filter response 
function as 

77(i,cj;r) = e(t)re-(r+^'^)* (17) 



where Q{t) is the unit step function. 

Figure[S]shows the real part of the two-time correlation 
function 

Git,t + T) ^ {ali^{t)a^k{t + T)). 

for d — X/4 and in the range 0<t + T<T — 2TT/g when 
the atoms are in a Mott insulator state, while Fig.[n]shows 
that same function for the initial superfiuid state IV'SFi), 
and Fig. [7] for the initial superfluid state |V'sf2)- As was 
the case for the reflected intensity, the periodic oscilla- 
tions of the reflected intensity characteristic of the Mott 
insulator case are washed out in the superfiuid regime. 
Since the number-conserving description |V'sf2) of the su- 
perfluid state corresponds to a smaller atom number un- 
certainty, the corresponding correlation function is char- 
acterized by stronger oscillations, intermediate between 
its description in terms of coherent states and the Mott 
insulator case, see Fig. [7] . 

In terms of r = ti — t2, the physical spectrum S{t, lo] V) 
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Re[G{t + T, t)] 




Figure 7: Two-time correlation function G{t + r, t) for atoms 
initially in the superfluid state \'4>sy2) with a total number of 
8 atoms, d — A/4 rj = O.lg and 7 = 0.5g. The curve is an 
over 50 trajectories at time t and 3 trajectories at time r. 



Re[So{io)] 




Figure 8: Real part of So{u}) for a Mott insulator states with 
6 atoms in well and 2 atoms in well 1 , with d — A/4. The 
peak is centered at \no — ni| in this case. 



may be reexpressed as 

Jo 

nt-r 

X / dtse^"^ Re[(a^_^(i2+T)a_fe(i2))] 
= 2r2e-2rt5'p(^^r). (18) 

with r > 0. 

Figure[5]shows the real part of Sq{ijj) for atoms initially 
in a Mott insulator state, while Figs. [9] and [10] are for 
atoms initially in superfluid states. In Fig. [S] the tem- 
poral oscillations of the correlation function G{t -I- t, t) 
translates into a modulation of the lorentzian spectrum 
due to the filter function at the oscillation frequency 
Lo^ = g\nQ — ni \ between the incident and reflected light 
fields. The physical spectrum provides therefore a di- 
rect measure of the population difference between the 
two wells for d = A/4, of their sum for d — A/2 and 
hence of no and ni separately, i.e., 



(19) 



where we assumed no > ni for concreteness. The sec- 
ond small dip at frequency g{no + rii) originates from 
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Figure 9: Real part of So{i^) of superfluid states IV'sfi) with 
each well containing 4 atoms on average and d — A/4. 




Figure 10: Real part of So{ui) for the superfluid state |?/'sf2) 
with a total number of 8 atoms and d = A/4. 



the pump term, /i77(a|,e'^"''"'""^-** -I- h.c), in the Hamilto- 
nian PT|) . 

Figures [9] and \W\ show the physical spectrum for the 
coherent and number-conserving superfluid states iV'SFi) 
and |V'SF2), respectively. In the number-conserving de- 
scription (I14[) the total atom number N is fixed, hence for 
d = A/4 the population difference between the two wells 
|2no — A^l can only take all even or all odd numbers, de- 
pending on N. This restriction disappears in the coher- 
ent state description of IV'SFi), so that additional peaks 
appear although the mean number of atoms in each well 
are the same in Figs. [HI andfTUl . 

We conclude this discussion by mentioning that for a 
well separation d = X/2 case the superfluid descriptions 
|?/'SFi) and \'ipsF2) result in completely different spectra. 
For IV'SFi) the spectrum is similar to that of the d = 
A/4 case since no and ni are independent and hence give 
all combinations of different total numbers (no + ni) to 
form the spectrum. By contrast, the number-conserving 
description |'(/'sf2) results in a single sharp peak, just as 
in the Mott insulator case, see Fig. [51 

Next we use the number-conserving approximation 
|V'SF2) to investigate the relation between the full width 
at half maximum(FWHM) of the physical spectrum and 
the atomic number uncertainties in the two wells. We 
evaluate that width with the help of an envelope func- 
tion fltted as illustrated in Fig. [TlJ The atomic number 
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Figure 12: FWHM w vs. atomic number uncertainty of each 
well a. The data points corresponds to A''=l,2,3,4,5,6,7,8,10 
and 12. The straight line is the fit to the data points. 



uncertainties a are defined in the usual way as 



\ 



N 

E 

m=0 



(to 



62 



(20) 



where to — X^m'^^m' ^ ^^'^ total number of atoms, 
and 6^, the probability of having to atoms in one of the 
wells, is a binomial distribution. Figure [T^l which plots 
the FWHM w of the physical spectrum as a function of 
tJ, shows as expected a linear dependence, further illus- 
trating the use of that spectrum in helping characterize 
the many-body state of the atomic system. 



V. ENTANGLEMENT 

In this section we make some brief remarks on the en- 
tanglement generated between the light fields and the 
atoms, and its dependence on the state on the atomic 
system. We characterize the entanglement in terms of 
the logarithmic negativity fl9j 



Em{p) = log2 Hp 



(21) 
1 denotes 



where p^^ is the partial transpose of p and 
the trace norm. 

There is always some ambiguity in the way a quantum 
system is described in terms of subsystems. One simple 



way to describe the system at hand is as a four-partite 
system comprised of two optical modes and the atoms 
in the two wells. It is then possible to trace over any 
two parts of the full system, and to consider the remain- 
ing subsystems only. We find that independently of our 
choice of the subsystems being considered the resulting 
reduced density operator does not retain any trace of 
whatever entanglement may have characterized the full 
system. This leads us to describe instead the full sys- 
tem as a bipartite system, the two subsystems being the 
optical field and the atoms. 

The situation is particularly simple if the atomic sys- 
tem is initially in a Mott insulator state. As discussed 
in Section [TTTl in the absence of quantum jumps between 
the times and t we have that 

\ip{t)) = |afc(no,"-i,c',i),a_fc(no,ni,d,t);no,ni) (22) 

where the amplitudes of the coherent states au and 
a-k depend on the atomic populations of the wells and 
on their separation. It is immediately apparent from 
Eq. that the total optical field is not entangled with 
the atoms. As previously remarked, in the case of co- 
herent states, which are eigenstates of the annihilation 
operator, quantum jumps do not affect the light field, 
hence \'>p{t)) is always a product state for this four-partite 
system, that is, no entanglement builds up between the 
atoms and the field. 

The situation is markedly different for atoms in a 
superfiuid state. For example, in the atom number- 
conserving approximation |V'sf2), and in the absence of 
quantum jumps between times and the state of the 
system at time t is 

N 

IV'(O) = '^bno\ak{no,d,t),a-k{nQ,d,t)\no,N - uq), 

no 

(23) 

a state for which the light field and atoms are clearly en- 
tangled. Dissipation and the associated quantum jumps 
will also clearly impact the amount of entanglement in 
that case. To illustrate that point, we first evaluate the 
logarithmic negativity of the bipartite atoms-optical field 
system for a well separation of d = A/4, while neglecting 
pump and decay mechanisms. In that case, from Eqs. ([3]) 
and (O the system is described by the Hamiltonian 

H = %[(no + "-i)(aiafe + aLfeft-fc) 

-I- (no - ni)(aLfcafc + ala-fc)], (24) 
and we find easily that 
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where 



(Xk{no, ni,t) 
a-k{no, ni,t) 



a , _ 

a , . 
2^^ 



■2ignot _|_ ^~2ignit\ 



-2ign(,t 



-2ignit\ 



(25) 
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Figure 13: The logarithmic negativity of light and atoms for 
a perfect cavity(no decay and no pump) . The initial coherent 
amplitude of k mode a — 1 and atoms are in superfluid states 
|V'SF2) with total number of atoms N — 8 and well separation 
d = A/4. 



As a final point, we remark that similar results hold 
in case the superfluid atomic system is described by the 
state IV'SFi), although the maximum value of the loga- 
rithmic negativity is now larger, due to the larger num- 
ber of terms in the expansion of that state in terms of 
number states. 
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Figure 14; Logarithmic negativity for a perfect cavity, with- 
out decay but including the pump. The initial coherent am- 
plitude of the k mode is a = 1 and the atoms are initially in 
the superfluid state |V'sf2) with total number of atoms N = 8 
and well separation d = A/4. The pump constant is r; = O.Sg. 



The logarithmic negativity Ej^ corresponding to that 
case is shown in Fig. [T3] for = 8 atoms. We observe 
that E}^ = for t — iir/g where € is a integer, indicat- 
ing the absence of entanglement between the light and 
the atoms. These are the times when the atoms do not 
reflect any light into the — fc-mode, and the state of the 
system has undergone a full revival to its original, unen- 
tangled form. 

The situation is changed when including the external 
pump in the description of the system, in which case there 
is no longer an exact revival. Figure [2] illustrates that 
a remnant of that feature is still observable in that case, 
but is disappears almost completely when both pump and 
dissipation are included, see Fig. [13 In that latter case, 
the entanglement between the light field and the atoms 
disappears completely over time as would be expected. 



Figure 15: Logarithmic negativity including dissipation. The 
initial coherent amplitude of the k mode is a = 1 and the 
atoms are in a superfluid state \'ipsF2) with total number of 
atoms N = 2 and well separation d = A/4. The pump con- 
stant is = Q.2g and the decay constant is 7 = 0.7g. The 
curve has been averaged over 100 trajectories. 



VI. CONCLUSIONS 

We have studied the interaction of a light field and ul- 
tracold atoms either in the Mott insulator state or in the 
superfluid state in a two-well optical lattice in a lossy 
cavity. We have calculated the reflected intensity, the 
physical spectrum and the logarithmic negativity, and 
shown that all three observables present completely dif- 
ferent features for the two different manybody atomic 
states, allowing one can distinguish these states. 

Future work will generalize this analysis to more 
complex manybody atomic states, using in particular 
Laguerre-Gaussian modes of the light field to character- 
ize the state of vortex lattices. We are also extending this 
work to optomechanical situations, using radiation pres- 
sure force on Fabry-Perot resonators with moving mirrors 
to induce quantum phase transitions and monitor them 
optically in real time. 
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